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Based on a recent progress in nonequilibrium statistical mechanics of infinitely extended 
quantum systems, a nonequlibrium steady state (NESS) is constructed for a single-level 
quantum dot interacting with two free reservoirs under less general but more practically 
useful conditions than the previous works. As an example, a model of an Ahoronov-Bohm 
ring with a quantum dot is studied in detail. Then, NESS is shown to be regarded as a 
MacLennan-Zubarev ensemble. A formal relation between response and correlation at NESS 
is derived as well. 

§1. Introduction 

Nonequilibrium statistical ensembles have been studied for many years, but no 
consensus has been made. As an illustration, let us consider a MacLennan-Zubarev 
ensemble for a system consisting of identical particles. The system is assumed to be 
divided into M independent parts, each of which has the energy Hj and the particle 
number Nj (j = 1, • • • M) and which are interacting by an interaction W. According 
to the MacLennan-Zubarev approach,^)' a steady state close to a local equilirium 
state is described by: 

P+ = ^ = -exp[-^f3jiHj-f,jN,)+ dse^'Jsis)}, (M) 

j=i 

where Z is the normalization constant, and fij are, respectively, the temper- 
ature and chemical potential of the jth. subsystem, Hj = Hj -\- dse*"^ '^^dl^^ ' 
and Nj = Nj -\- dse'''^ '^^df'^ ^'^^ Zubarev's local integrals of motion,^) and 

Hj[s) = e^Hs/hjj.^-iHs/h^ ^ ^iHs/hjy.^-iHs/h ^^^j^ = ff. _^ ^ ^j^g ^^^^^i 

Hamiltonian. The integrand in the left-hand side is given by Js{s) = "^fLiPjJjis) 

where Jj{s) = '^^^^^ — l^j~~i^ stands for non-systematic energy flow, or heat flow, 
to the jth subsystem. A convergence factor e''^ (e > 0) is introduced in the time 
integral, where the limit e — > is taken after all the calculations. As discussed 
in Ref. 2), this ensemble well describes nonequilibrium phenomena, but it has a 
fundamental difficulty. Indeed, because Js{s) is the sum of heat flows divided by 
subsystem tempertures, it is the entropy production rate of the whole system. Hence, 
if the ensemble would describe a state consistent with the second law of ther- 
modynamics, the average of Jsis) over p+ should be a positive constant and, as a 
cosequence, the integral in (jl-lj) would diverge in the limit of e ^ 0. 



typeset using VTpT^.ds (Ver.0.9) 



2 



S. Tasaki, J. Takahashi 



On the other hand, rigorous researches have been carried out on nonequihb- 
rium steady states (NESS) of infinitely extended systems and are developed further 
in recent years. Those include studies on NESSs of harmonic crystals, ^^'^^ a one- 
dimcnsional gas,^) unharmonic chains, ^-^ an isotropic XY-chain,'')'^) systems with 
asymptotic abelianness,^) a one-dimensional quantum conductor,^'') an interacting 
fermion-spin system,^^) fermionic junction systems, a quasi-spin model of super- 
conductors,^^^ a Bose-Einstein condensate in a junction system, a Bose-Einstein 
condenstate in a small system coupled with a large reservoir, ^^-^ a quantum dot cou- 
pled with several reservoirs, on nonequilibrium entropy productions^^-*"^'^^ and on 
linear responses.^)' See also reviews. ^^^"^^^ Moreover, we have shown^^^'^^^ that 
NESS constructed by this method well explains experiments on transports of some 
mcsoscopic systems. 

In this article, we illustrate the above-mentioned features in terms of a spinless- 
electron model of a single-level quantum dot interacting with two two-dimensional 
reservoirs and show that NESS can be regarded as a MacLennan-Zubarcv ensemble 
in an appropriate sense. In the next section, the basic features of the C*-algebraic 
method are summarized in a less technical way. In Sec. 3, a nonequlibrium steady 
state (NESS) is contsructed starting from a local equilibrium state with the aid of 
the scattering theory. Since general results were proved by Ruelle^-* and by Frohlich, 
Merkli, Ueltschi,^^) here we describe the construction under restricted but practically 
useful conditions and apply it to a model of an Ahoronov-Bohm ring with a quantum 
dot.25) In Sec. 4, wc show that NESS is an analog to the Zubarev-MacLennan 
ensemble under slightly different conditions from Rcf. 23). As an application of this 
observation, we study a formal relation between response and fluctuation at NESS. 
The last section is devoted to concluding remarks. 

Before closing this section, we describe the spinless-electron model of a single- 
level quantum dot interacting with two reservoirs. The system is described by cre- 
ation (annihilation) operators of the reservoir electron with wave number k G R^: 
a|.j.(ajkr) (r = L: for the left reservoir, r = R: for the right reservoir) and the dot 
electron: c^{c), which satisfy the canonical anticommutation relations: 

[akr, «Iv]+ = Srr'S{k - k')l , [c, c^]+ = 1 , 

where [A, B]^ = AB + BA is the anticommutator, Srr' the Kronecker delta, S{k — k') 
the Dirac delta function, 1 the unit operator and the other anticommutators vanish. 
The total Hamiltonian is H = Hl + Hr + Hd + W where 

Hr = j dkujkral^akr (r = L or R) (1-2) 

are the reservoir Hamiltonians, Hd = eoc^c the dot energy and W stands for the 
interaction among the reservoirs and the dot. In the above, ojkr {r = L, R) stands for 
the single-electron energy of reservoir electrons: lo^l = ^kR — ^y = j {2m) — eV/2 
where V is the bias voltage, e the elementary charge, h the Planck constant and m 
the effective mass, and eo is the energy of the dot level. The numbers of particles in 
the reservoirs and the dot are given, respectively, by iV^/i^ = / dA;a|.^y^afc£,/jj and 

Nd = c^c. 
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§2. C*-algebraic approach 

2.1. C* -algebra and Time Evolution^'^^'^^^ 

An essential feature of the C*-algebraic method is to discuss the properties of 
infinitely extended systems through the investigation of finite observables. We start 
from a set J- of operators A such that the maximum eigenvalue (more precisely, the 
maximum spectrum) of A^A is finite and its square root, denoted as ||^||, is used 
for measuring the size of A (or || • || is a norm). The set J- is a complex linear space 
where the product and the 'conjugation' A ^ A^ are defined*^, and the norm || • || 
satisfies (i) ||^|| > and ||^|| = implies ^ = 0, (n) \\aA + B\\ < \a\\\A\\ + ||_B|| 
{aeC,A,B e F), (iii) \\AB\\ < \\A\\\\B\\ and (iv) the C*-property: \\A^A\\ = \\Af. 
Also J- is complete with respect to this norm**-*. Such is called a C*-algebra.^'^^'^^^ 

For the spinless electron model of a quantum dot, ^ is a set of operators which 
can be approximated, with arbitrary precision, by a finite sum***-*: 

al + J]C7c6i(C)^---6,(C)Vi(C)---H(C) , (2-1) 
C 

where a, are complex numbers, 1 G T is the unit and bj{() (j = 1, • • • A^^) is either 
c or ar{f) = J dkf{k)*akr {r = L,R) with f{k) a square integrable function over 
R^, i.e., / G L^(R^). Namely, the algebra is generated by l,c,aL{f) and aji{f). 
Because of the canonical anticommutation relations, one has [a,.(/), 0^/(5)]+ = 0, 
[ar{f),ar'{g)^-i- = {f,g)6rr'^ where {f,g) = f dkf{k)*g{k). The definition of T is 
meaningful since ar{f) is bounded. Indeed, the C*-property leads to 

l|ar(/)r = ||a,(/)ta,(/)|p = ||a,(/)ta,(/)a,(/)ta,(/)|| 
= \\ar{f)H{fJ) - a.(/)ta.(/)}a.(/)|| = (/, /)||a.(/)f 

or ||ar(/)|| = \/ if, f) < +00. For later use, we introduce a subset .7>es(C of 
reservoir operators, each element of which is approximated, with arbitrary precision, 
by a finite sum: al+Ec C'c"n(/c,i)^ ' ' ' ar,(/c,s)^ar,+i(/c,s+i) • • • ^r^v^ ifcN^) i^j = L 
or R, /^ j G L^(R^), and a, are complex numbers). 

Now we turn to the description of the time evolution. Because of their un- 
boundedness, the Hamiltonians i?L, Hr are not included in T and, thus, are not 
observables. But, the time evolution can be defined within the framework of J-. 
As an example, let us consider the free evolution g^^ot/fi-^^-iHot/h gg^erated by 
Hq = Hl+Hr+Hd. a formal calculation gives e^"°^/^ar{f)e-^"°^/^ = f dkf{k, tfakr 
where f{k,t) = f{k)e'^^-^/^ G L2(r2), and e^"ot/h^^-iHot/h ^ e-'^^^/^c, both of 
which are elements of J-. Thus, one can define the time evolution of the finite sum 
(|2-1|) and, thus, of any element of J^. Hereafter, to avoid the explicit use of Hq, the 

*' In the mathematical literatures such as Ref. 27), A* is used instead of A'' . 
**^ Namely, if {A„}55!Li C .7-" is a sequence such that \\An — Am\\ for n,m 00, then there 
is A £ such that \\An — A\\ — > as n ^ 00, or {j4„}5JLi has a limit in JT. 

***-* Let A G T, then, for arbitary e > 0, there exist a, C( and tij(C) {j = 1, ■ ■ '-^c) such that 
\\A - {al + CMCV ■ ■ ■ bsiCVbs+iiO • • • K (0}ll < E holds. 
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evolution is denoted as rj; ' (A) = e^Hot/hj^^^iHot/h^ ^^^^ ^^le map A rf' (A) is lin- 
ear and preserves the product, conjugation and norm: Tj:^\AB) = Tj:^\A)T^^\B) , 
T^^\Ay = t^^\a'') and ||rj''°'*(^)|| = ||^||. Its infinitesimal generator is given by 
d{A) = 4j:T^'^\A) on some dense subset D{6) C called the domain of 6*\ 



t=o 



For instance, 6(ar{f)) = —i/h j dkujkrf{k)*akr (r = L,R) is meaningful only when 
^krf{k) € L^(R^), and a set of finite sums generated by such ar{f) (r = L,R) to- 
gether with 1 and c provide the domain D{5) of 5. The evolution Tt generated by 
the total Hamiltonian H = Hq + W is defined as a solution of 

|r*(A) = Tt (^6{A) + ^[W, A]^ , CA e D{S)) , 

under the initial condition Tt{A)\f=Q = A. The map Tf has similar properties as r^^-*. 
2.2. States^V,2S) 

Usually, a statistical state is given by a density matrix. However, within the 
algebraic approach, it is specified by listing the average value {A) of an arbitrary 
element A £ J^, i.e., by a complex- valued linear map, called a normalized positive 
linear functional: A {A), which satisfies (A'^A) > 0, (1) = 1 and \ {A)\ < \\A\\. 

Within the algebraic approach, canonical states are formulated without explicit 
reference to the Hamiltonian. Remind that the grand canonical state with temper- 
ature and chemical potential jU of a finite degree-of-freedom system is given 
by {A)gc = Ti:{Ae-^'^^->^^'>)/Zgc with H the Hamiltonian, N the total number 
of particles and Zgc the normalization constant. Then, it is easy to see that the 
Kubo-Martin-Schwinger (KMS) boundary condition^o)'^!) {Aafl^{B))gc = {BA)gc is 
satisfied with respect to as'^{A) = QAH-tiN)s ^^-i{h-ij,n)s _ Yqi infinite systems, 
the KMS condition defines canonical ensembles. ^^^'^^^ As an example, let us con- 
sider the grand canonical state with temperature and chemical potential jU of 
the reservoir system described by J^rcs- Let ai^{ar{f )) = / dkf{k)*e~'^^^'^''~^^^^akr 
(r = L, R) and J^iei'^ be a dense set such that, for any A G J^ei'^, ai'^{A) is analytic 
in |Ims| < then the grand canonical state is defined as a state satisfying 

{Aaf^{B))g, = {BA)g, , {A, Be J^^i^) . 

This equation and the canonical anticommutation relation lead to 

{ol{f){'^i'^{o-r'{9))+ar'{9)})gc = (KKf), 4(/)]+)gc = Srr' {q, f) = Srr' J dkg{k)* f (k) 

Since afp{ar>{g)) + ar'{g) = f dkg{ky{el^^^kr'-t^) + 1)0^^,, by replacing g{k) by 
g(k)F{ukr) with F{x) = l/{e^(^~'*) -1-1} the Fermi distribution function, one obtains 

{al{f)ar'{g))gc = 5rr' j dkg{k)*f{k)F{ujkr) ■ 



*•* Namely, any A £ can be approximated by an element of D{6) with arbitrary precision. 
**) .T^e'r can be a set of finite sums generated by o,.(/) with f{k)e''^'^>"-i^^'' eL'^(R^) (|Im s\ < P). 
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In the same way, one can show that the state {• • •)gc satisfies Wick's theorem. In 
short, the KMS condition fully determines the state {■ ■ ■)gc- In general, if no phase 
transition takes place, the KMS state is unique and, if a phase transition occurs, 
several KMS states exist as a result of spontaneous symmetry breaking. 

A local equilibrium state can also be defined as a KMS state. As an example, we 
consider a local equilibrium state where the left (right) reservoir is in the equilibrium 
state with temperature (3]^^ (Pr^) chemical potential m (/U_r)- Consider a 
map as formally expressed as as{A) = e^T.rMHr-iirNr)s^^-iY.r3r(Hr-tirNr)s 
defined by (7,(a^(/)) = J dkfik)*e-'^'~^'^kr-i^r)s^^^ ^ and (t,(c) = c, then a 

local equilibrium state (• • • );oc is given by a KMS condition {Aai{B))ioc = {BA)ioc 
{A e T,B e J^res)' where B G T^^^ implies that (Js{B) is analytic in |Ims| < 1. It 
again satisfies Wick's theorem and its nonvanishing two-point functions are 

{al{f)ar{g))ioc = Jdkg{kyfik)Fr{uJkr)ir = L,R), (c^^oc = ^ (2-2) 

where Fr{x) = l/{el^rix-f^r) + 1} (r = L,R) are the Fermi distribution functions. 

2.3. Ergodicity^'^)'^^) 

If the decay of d3mamical correlations is sufficiently fast, certain states have 
ergodicity. One of such dynamical conditions is the asymptotic abelian property:^^) 

lim \\[A,rt{B)]\\ = , {A or B e J'even) (2-3) 
lim \\[A,TtiB)]+\\ =0, {AandBG ^odd) (2-4) 

|t|— +00 

where ^even/odd = {A E : A consists of even/odd numbers of Fermion operators.} 
and [A, B] stands for the commutator: [^4, B] = AB — BA. Then, as stated below 
(Example 4.3.24 of Ref. 27)), a clustering property is satisfied by a class of states 
called 'factor' states, which include unique KMS states. Roughly speaking, a state 
(• • • ) is called 'factor'*) if any D e T satisfying {[A, D]) = (VA G T) behaves as 
some complex number 7 in the sense of (ADB) = ^{AB) {^A, B E J^). 

Clustering Property:^^) For an asymptotic abelian evolution Tt and a 
factor state (•••), one has 

hm {{ATt{B)C)-{AC){Tt{B))} = 0. 

\t\^oo 

If {Tt{A)) = {A) (^A e J^), it is mixing: lim\t\^^{Art{B)C) = {AC){B). 

For the spinless electron model of a quantum dot, we show that the local equilibrium 
state {■■■)ioc restricted to the subalgebra J^j-cs generated by \,ai{f) and aR{f) is 

*^ Given a state (• • • ) over a C*-algebra J^, it can be represented as a subalgebra n{J^) of 
the algebra B of all bounded operators on some Hilbert space Ti. as {A) = {Q,-K{A)fi) with a 
cyclic vector f2 £ Ti (GNS representation). Let = {a £ B : [a,b] = 0,V& € 7r(-^)} and 

tt{J^)" = {ceB:[c,a] = 0, Va £ 7r(J'^)'}, then (• • • ) is called a factor state iff 7r(.F)' n n{J^)" = CI, 
where C is the set of complex numbers and 1 is the unit of 7r(.F). 
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mixing with respect to the evolution \ The stcite (* * * );oc 

is a factor state as 

a unique KMS state and is easily shown to be r^'''^'* -invariant. On the other hand, 
for A,B£ J-res, their commutator [A,t^^\b)] can be approximated with arbitrary 
precision by a finite sum of the terms like Ct[ar{f),T^^\ar{g)^)]+Dt {Ct,Dt G .^res) 
and its conjugate. Then, the Riemann-Lebesgue theorem^^^ gives 



Ct[arif),T^'\ar{g)^)]+Dt\\ < sup[||Ct|||| A| 



oo) 



Thus, one has \\m.\t\^^ \\[A,tI^\b)]\\ = 0, or rf^"* is asymptotic abelian. Therefore, 
as a result of Clustering Property, local equilibrium state (• • • )/oc restricted to .7>es 
is mixing with respect to t^^^: 

hm (^rf ) {B)C)ioc = {AC)ioc{B)ioc ■ {A, B,C e .F^es) (2-5) 

>oo 

Note that the state (• • • )ioc is not ergodic on the whole algebra with respect to 
T^^^ because {c^ Tj:^\c)) loc = e.~^^°^/^{c^ c)ioc does not converge for \t\ oo. 

§3. Nonequilibrium Steady States 

3.1. Scattering Problem and Nonequilibrium Steady States 

As discussed in Refs. 3), 6), 7), 9)-14), 18), 19), 22)-24), a nonequilibrium steady 
state (NESS) (• • • )± is constructed dynamically as an asymptotic state starting from 
the local equilibrium state (• • • )/oc introduced in Sec. 12.21 {A)± = limt^±Qo{Tt{A))ioc- 
As discussed by Ruelle^^ for systems with L^-asymptotic abelian property and by 
Frohlich, Merkli and Ueltschi"'^'^) for quantum junction systems, the construction of 
NESS is closely related to the scattering problem. Here we give a restrictive but 
practically useful characterization used in Refs. 25), 26). 

For the spinless model of a quantum dot, the interaction W induces scattering of 
the left /right-reservoir electrons and the process is described by asymptotic fields:^*^^ 



ar/°"^^(/)= lim Tt[Tr-'{ar{f)) , {r = L,R) (3-1) 

where °^°"*^(/) are incoming/outgoing fields and the limit is taken in an approriate 
sense. Here we consider a case where the limit ()3-l() exists in norm: 



lim 

— 00/+00 



Tt(rf)-\a,(/)))-a(-/-*)(/) =0, {r = L,R) (3-2 



and the initial state (• • • )/oc is r^'''^'* -invariant. Then, one has^^^'^^^ 

Proposition 1: If, for some subset Dq C L^(R^), (i) the limits (|3-2() exist 
for any / G Dq, (ii) the fields a^"^^°''^\f) {r = L,R, f e Dq C L2(r2)) 
generate the whole algebra JF and (iii) (• • • )/oc is r^'^^ -invariant, then, the 



limit 



hm {n{A))ioc, (Me.F) (3-3) 
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exists and defines a (nonequilibrium) Tf-invariant state (•••)+/-• More- 



over, 



{aiT\h)^ • • • a(r')(/.)ta(°f^)(/,+i) • • • ^r^C/^))- 
= (a(^)(/i)t • • • a^^^\fs)^ai%{U^) ■ ■ ■ a^;^\Mh 
= (ari(/i)^---ar,(/s)^ar,+i(/s+i)---arjv(/Af))ioc • (3-4) 

N.B. For a model of a single- level quantum dot coupled with free reservoirs (SEBB 
model) where the interaction is bilinear with respect to field operators, Aschbacher, 
Jaksic, Pautrat and Pillet^^^ derive an equivalent characterization to Proposition 1, 
but the present form is applicable, in principle, even to the case where the intraction 
is not bilinear. 

In the previous section, we have seen that the unperturbed evolution restricted to 
the reservoir algebra J-^es is mixing. As a consequence, as first shown by Ruelle^^ (see 
also Ref.24)), the steady states (• • • )± mixing with respect to the full evolution 
Tt in both directions of time: 

Proposition 2: Under the same conditions as Proposition 1, the steady 
state (• ••)+/- satisfies 

lim {ATtiB)C)+/_ = {AC)+/_{B)+/_ ,iA,B,CeJ^). (3-5) 

Proof of Proposition 1: It can be shown immediately from the following lemma: 

Lemma 3: If the limits (|3-2() exist and the fields ar^^°"^\f) (r = L,i?) 
generate the whole algebra there exists ^j^i_{A) G T for any A ^ T 

such that limt^_oo/4.oo \\^ — 'Tt{Tf^'^ ^(7+/_(j4))) || = 0. The map 7± is a *- 
isomorphism between T and the reservoir algebra J-^es generated by ar(/) 
(r = L,R), namely, (i) 7±(A) € J^^es, (h) l±{aA + B) = a7±(^) + 7±(^), 
(hi) 7±(j4i?) = 7±(j4)7±(i?), (iv) 7±(j4^) = '^±{A)'^ and (v) 7± is one-to- 
one and onto. Moreover, (vi) \\A\\ = ||7±(^)|| and 



7- • • • 4r')(/s)^4-*)(/,+i) • • • a(7)(/^)^ 

= 7+ (<^(/i)^ • • • a^:\f,)^a^fl{fs+i) ■ ■ ■ 4;^{fN\ 
= ari(/i)'^---ar,(/s)^ar,+i(/s+i)---arjv(/Af) • (3-6) 



The maps 7± are nothing but the M0ller morphisms:^-*''^'^^''^'^^'^^^ 

hm ||rf-^(r,(^))-7+/-(^)||=0. (3-7) 

t— »+oo/— oo 



Indeed, this lemma and Tj^'^^-invariance of {■ ■ ■)ioc give 
{rt{A))ioc-{l±{A))ioc\ = \{n[A - r^t{rT\7±{m)) 



loc 



< 
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n[A - r_,(ri^\7±(^)))) || = P " r_,(ri^\7±(^))) II - (t - ±00) 



or {A)± = limt^±oo {n{A))ioc = (7±(^))«oc exists. Clearly, the map A {A)± is 
linear and, as {A^A)± = lim {Tt{A)^ Tt{A)) > and (1)± = lim {rt{l))ioc = 1, 

I— >±oo t— > + oo 

{• • •)± defines a state. And it is invariant: {Tt{A))± = limti-^±oo{Tt+t' (^))/oc = (^)± • 
By substituting A = a^^/'''''\h)^ • • • ai':/'''''\fs)U^Tir\fs+i) ■ ■ ■ aiT'^'hfN) into 
{A)± = {'y±{A))ioc, (EH immediately follows from (ESI). {Q.E.D.) 

Proof of Proposition 2: We only show it in case of {■ ■ Eq. (|3-7|) gives 

||rf (7+(^)) -7+(rt(^))|| < ||rf (7+(^)) - ri°) (rW-i(r,(^))) || 

+ \\rf};\T,^t[rt{A))) - ^^{n{A))\\ ^ - +00) , 

i.e., Tf^\'^+{A)) = 7+(rf(j4)) (^t G R). On the other hand, as shown in the proof of 
Proposition 1, we have = {'^j^{A))ioc and 7+(vl) G ^res-, thus, l|2-5p gives 

hm {ATt{B)C)+ = hm {^+{A)^+{Tt{B))^+{C))i,, 

\t\^oo \t\~*oo 

= hm (7+(^)Tf (7+(i?))7+(C))io, = {^+{A)^+{C))Ul+{B))ioc 
= {AC) + {B)+ . {Q.E.D.) 

Proof of Lemma 3: Let us consider the case of incoming fields. Because of ()3-2|) and 

\\Tt{Tt^^~^{ar{f)^)\\ = ||ar(/)"''||, One has 

\\n{ri''>-\ar{f )^aA9)))-a^rHf)^4^\9)\\ < ||a.(/)t|| ||r,(rf)-\a,K5))) " 4^^(5)11 

+ ||Tt(rf )-\a,(/)t)) - 4-)(/)t||||4r)(5)|| ^ ^ -00) 

Repeating the same arguments, one finds 

lim ||ri(ri°)-\a,,(/i)^^ • • • a,^(/^)^-)) - 4^\f^)^^ ■ ■ ■ 4^\fr,)^'^\\ = , (3-8) 

where \]j (j = 1, • • • A^) stands for f or no symbol. Therefore, for any finite sum: 

Af^al + Y. C(o!^^ (/c.i)^ • • • ifcs+i) ■ ■ ■ <J (ZcTVc) , (3-9) 

C 

where a,C( G C, there exists 

= "1 + H CC^nl/ci)^ • • • arAf(,s)^ar^+iif(,s+l) ■ ■ ■ arjv^ if(,Nc) e -^res (3-10) 

c 

such that limt_>_oo ||^/ — Tt{jt^^ ^(-^/))ll ~ 0. On the other hand, as the fields 
Or™^(/) {f = L,R) generate IF, any A G T can be approximated by a finite sum 
(|3-9() with arbitrary precision. Thus, for any A G F, there exists B G J-^cs such 
that lim ||j4 — tAt^^^ ^(i?))|| = 0. As such B is unique, we define B = '^^{A). 

t^ — cxs 
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Then, the properties (ii), (iii), (iv) and (vi) can be shown immediately. For example, 

(vi) follows from ||7+(^)|| - \\A\\ < \\Tt{T^^^'\j+{A)) - A\\ ^ {t ^ -oo). The 

property (vi) implies that 7+ is one-to-one. Moreover, since any B G J-^cs can be 
approximated by a finite sum H3-1U() . one can find A £ such that J+{A) = B or 
7+ is onto. The second equality of follows from the definition of 7+ and (|3-8j) . 
The properties of 7_ can be proved in the same way. {Q.E.D.) 

3.2. Nonequilibrium Steady States for Ahoronov-Bohm Ring with Quantum Dot 

In this subsection, we further investigate the properties of NESS (•••)+ when 
the reservoir-dot interaction is described by a sum of two tunneling interactions^^) 

W = J dk{ukLa\^c + UkRa\^c} + -we*"^ j dkdq UkLUqR o^kh'^qR + (h-c.) , (3-11) 

where the first term corresponds to a tunneling via a quantum dot and the second 
to a direct tunneling between the two reservoirs. Real parameters w and (/? are, 
respectively, the relative strength and phase between the two processes. This model 
describes an Ahoronov-Bohm (AB) ring with a quantum dot^^^ and, when t/; = 0, it 
reduces to a model of a single-level quantum dot embedded between two reservoirs 
studied in Ref. 24). The tunneling matrix elements are assumed to satisfy: 

(a) The real- valued functions Ukr (r = L,R) are infinitely differentiable with respect 
to A; G and Ukr = when \k\ < ko or \k\ > ki (for some < /cq < ki). Then, 
the functions rr{oj) = 27r J dk\ukr\^5{ijj — Wkr) {r = L,R) are integrable and 
infinitely differentiable on the whole real axis. 

(b) Let Mr{z) = J dk\ukr\'^ / {z — tokr), then the function 

A{z) = (1 - w^Ml{z)Mii{z)){z - eo) - ^ Mr{z) - 2w cos (pMl{z)Mr{z) 

r=L,R 

has no real zeros and, hence, 1/ A±{uj) = lim 1/ A{lo it ie) is bounded. 

e>0,e^O 

3.2.1. Construction of NESS 

To derive incoming fields, it is enough to study the evolution T^tTt^\o'ro{f)) 
where f{k) is in a set Cg°(R^) of infinitely differentiable localized functions*^ as 
Cg°(R^) is dense in L^(R^). Then, thanks to the bilinearity of W, T-t{aro{f)) is 
written as 

r-t{aro{f)) = ^ ar'{ipr'{t)) +ilJc{tTc . 

r'=L,R 

From the equation of motion j-^T-t{aro{f)) = -'5(-^-i(«ro(/))) - '^-t(aro (/))], 
the functions ipr'{k',t) and •i/'c(i) are found to satisfy 

d 



ih—tpr'ik; t) = LOkr'tpr'ik; t) + Ur' (k) ^ipdt) + we^'^''' J dk'u-,{k')Tp-,{k'; t)|, 
ih^Mt) = eoMi)+ Yl j dk'ur^{k')^r'{k'-,t) , (3-12) 

r'=L,R-^ 



*•* More precisely, functions with compact support. 
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with an initial condition: ipr'ik^O) = 6r'rof{k), V'c(O) = 0, where L = R,R = L, 
ipL = —fR = ^- The hnear equations p-l'if) can be solved easily and one has 



A^{L0i^,-,){L0k'r' — ^^kr' + ^0) Aj^{ljOj^,-,){ujy-, — Wfcr' + ^0) j 

V'c(i) = dk — — — Fkr'it) (3-13) 

r'=L,R'' +\^kr') 

where A±{u}) is introduced in (b), i^{x) = 1 + M^{x){2w cosip + tiP'ix — eo)), 
Kr{x) = l+it;e*'^'-(x-eo), xf{x) = l+we'-'^'' M^{x) with M^{x) = lim M^(x±ie), 

e>0,e^O 



Fkr{t)=e 



^ (/t)-f— [dk'l ^^ (^kr)Urik')lpr{k') ^ Kr{uJkr)Ur{k')ll)f{k') 



with ^/'r(fc) = Vro/(^)- When /(fc) G Cq°(R^) and the conditions (a), (b) are 
satisfied, the limits like / dk h{k) / {x — uj^r ^ i^) = 1™ / dkh{k)/{x — ujkr^ie) 

J e-*0,e>Oj 

converge pointwise and uniformly with respect to x. As a consequence, T-tTt^^ («ro (/)) 
is again a linear combination 



r'=L,R r'=L,R 



where '4^^™}i{k; f), iproc{f) are derived from (|3-13|) by replacing Fkr{t) to Srrof{k)'- 



(3-14) 



, nr'o(fc)^ro(fc0*^4(a;fc'ro) 
^+('^fc'ro)(<^fc'ro - <^fcro + ^0) ' 

and Aipr'{k;t), Aipdt) are obtained from H3-13() by replacing Fkr{t) to 

Ahrit) = i^,rj.r.{.kr) I dk' ^^0(^0/(^0 .K.,-...)V. ^ 

with ??roro(2;) = ifnix)-, Vforoi^) — i^rnix)- Reminding lim e~*^*/ {x + iO) = —2'Ki5{x) 
in the sense of distribution, one expects AFkrit) ^ (t ^ +co) and, thus, 
r_trf (a.„(/))- aA^^f^Mf)) + ^l^^HfTc ^ 4^\f) , (t - +00). (3-15) 

r'=L,R 

Indeed, if the conditions (a) and (b) are satisfied and / G Cj^(R^), an argument 
similar to those of Refs. 14), 29) gives 

\\r-M'\arM)) - 4^Hm < E \\aA^A'm + \^Mt)\ 

r'=L,R 
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= Yl J J dk\Ai;r'{k,t)f + \A'iljc{t)\^ ,{t^ +00) . 

r'=L,R V 

Moreover, under the condition (a), the original operators are expressed as 

arif) = ^ a^r^\^r'r{f)) , = X] "'r^\v^r'c), (3-16) 

r'=L,R r'=L,R 

where the functions ip^'rik; f) and ipr'dk) are given by 



This imphes that the incoming fields generate the whole algebra J-. Then, because 
of Proposition 1 and the fact that Wick's theorem holds for the expectation value of 
a product of a.f (/) and ar'{f')'^ with respect to (• • • )ioc, we find: 

Proposition 4: For the model where the interaction is given by (|3-1H) . 
if the tunneling matrix elements satisfy the conditions (a) and (b), then 
the limit: lim(_^+oo('^t (^))/oc = (^)+ exists for any A £ and defines a 
steady state. Moreover, the steady state (•••)+ satisfies Wick's theorem 
with respect to the incoming fields ai'^V) introduced in (nTTll) and (1^^ . 
and the nonvanishing two-point functions are given by 



(a(';^)(/i)taM(/2))+ = {al{fi)ar,{f2))ioc = j dkh{kr h{k)Fr,{uJUr,) , 

(3-18) 

where Fr{x) = l/{e^''*^^' ^''^ + 1} (r = L,R) is the Fermi distribution 
function of inverse temperature (3r and chemical potential fir- 

A simple interpretation could be given to this result. Suppose that there exists an 
invariant vacuum state |vac.), then, due to the difference between the Schrodinger 
and Heisenberg pictures, the vector r_f (ar™''(/)^) |vac.) describes a one-particle state 
at time t starting from an initial state: ar™^(/)^|vac.). When i <C 0, one may 
regard r_( (ar°^(/)^)vac.) ~ t^^^ {ar{f)'^)\yac.) as a consequence of a relation: 

limi^_oo||ri°Vr(/)^) -r„t(a^^(/)^)ll = (cf. eq.(E3). In this sense, a^)(/)t 
describes a particle which was an unperturbed particle of the rth reservoir in the far 
past. Thus, (•••)+ is a steady state such that particles carry the temperature and 
chemical potential of the reservoir from which they come. 

3.2.2. Transports 

As an application of Proposition 4, transports in the steady state (•••)+ will be 
studied. Formally the energy and the particle number of the reservoir are expressed, 
respectively, hy Hr = J dkijOkra\^akr and Nr = f dka\^akr {r = L,R), and a formal 
calculation leads to 

ofJilh^^ ^- (af.(Ur)^{c + We*'^'' Of (Uf)} — (h.c.)) = Jr , 



-e'^^/^Nre-'^'/^ 
dt 



t=o h ' 
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^e'^'I^Hre-'^'/^ _ =-Uar{u^)^{c + we'^^af{ur)]-{h.c.)) =jf ,(3-19) 

where Ur{k.) is the tunnehng matrix elements, uf{k) = uJkrUr{k), and Jr,J^ G J- 
are defined by the middle expressions. Therefore, Jr and (r = L, R) can be 
regarded as the particle and energy flows to the reservoirs. Then, (|3-16jl . (|3-17|) and 
Proposition 4 give 

Proposition 5: The steady state considered in Proposition 4 carries the 
particle and energy flows: 

(Jl)+ = -{Jr)+ = ^I^ du;Tiu;){FR{u;) - , (3-20) 

1 f'°° 

(Jf )+ = -(J|)+ = —J ^du;T{o,){FR{u;) - F^M} . (3-21) 



where 

T{UJ) = -r—^^ rL{uj)rR{uj) . 

The entropy production rate Js = J2rl^rJr, where J? = — firJr are 
heat flows to the reservoirs, has the following NESS average 

{Js)+ = P duj^{FR{uj) - Fl{u)}{Pl{u^ - fiL) - Pr{uj - /iij)} , 

which is nonnegative and vanishes if and only if Pl = Pr and fj-i = fJ-R- 
The nonnegativity of {Js)+ immediately follows from an inequality: 

1 1 



ev + 1 + 1 



ix-y)>0 



where the equality holds if and only if x = y. As shown in Refs. 16)-19), positivity of 
the entropy production can be proved for more general cases since it is related to the 
relative entropy between the initial state and the state at time t: If the two states 
were described by density matrices, respectively, pioc and pt, the relative entropy 
S{pioc\pt) = Tr{/)((log/9t — log pioc)} is related to the entropy production via 

{Js)t = J^S{pioc\pt) , (3-22) 

where {■ ■ ■)t stands for the average with respect to pt- The same relation holds for 
a C*-generalization of the relative entropy given by Araki.^"^)'^^)"^^-* Then, because 
of S{pioc\pt) > 0, I'Hospital's rule gives the positivity of {Js)+^^^ 

For states described by density matrices, H3-22p can be easily shown. Then, 
we have pioc = e" ^^^^^"^^^^ V^o, Pt = e'^^^/^pioce^"^/'^ with Zq a constant, and 

S{piM = Tv{pi,,{ -Y^PoiHj - p,Nj) + Y,PjiH,{t) - PjN,{t)))} , 
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where Hj{t) = e^^t/^Hje"'^^/^ and Nj{t) = e'^^/^Hje''^^/^, and which gives (l3^ : 

j^S{pM = Tr{p^ocE/5.(^ -Z^.^)} = ^Js{t))ioc = {Js)t . 

j 

This observation indicates that the features consistent with thermodynamics come 
from the second term of the relative entropy: —Tr(/3j log p^oc)) which is similar to 
the nonequlibrium entropy of Zubarev:^) Sz = — Tr(ptlogpi) (cf. eq.(22.31) of 
Ref. 2)) with pi a reference local equilibrium state. An entropy of Zubarev type (more 
precisely, the relative entropy between the state pt and a reference local equilibrium 
state pi: S{pt\pi)) was studied by Frohlich, Merkli, Schwarz and Ueltschi^'^) in a 
slightly different context. 

The identification of Jr = ~ PrJr with the heat flow could be justified since 
dtJr{t) behaves as a thermodynamic heat in the weak coupling limit for a small 
system coupled with a single reservoir. 

Note that the expressions of the particle and energy flows agree with those of 
the Landauer formula. ^^^''^^^ As discussed e.g., by Sivan and Imry,^''^ when the tem- 
perature difference jSj^^ — {3^ and chemical potential differece pl — PR are small, 
(|3-20|) and (|3-21|) reduce to linear relations between thermodynamic forces and flows, 
where Onsagar's reciprocal relations hold. The general proofs of linear response rela- 
tions for nonlocal perturbations such as the temperature difference and/or chemical 
potential difference are discussed by Jaksic, Ogata and Fillet. 

§4. KMS Characterization of Nonequilibrium Steady States 

As discussed in Sec. 12.21 canonical and grand canonical states are character- 
ized as KMS states. As shown for L^-asymptotic abelian systems in Ref. 23), the 
nonequilibrium steady states discussed in the previous section can be characterized 
in a similar way: 

Proposition 6: If the limits (|S]3) exist, the fields a^r'''°''^\f ) (r = L, R) 
generate the whole algebra !F and {• • •)ioc is r^'''^^ -invariant, the nonequi- 
librium steady states {■■■)± are KMS states with respect to the maps 

aW=7±^a,7±, (4-1) 

namely, {Aa^{^\B))± = {BA)± for A, B G J^^ where i = (Tx IS a map 

defining (• • • )ioc as a KMS state, 7± are maps introduced in Lemma 3 and 
C J- are (dense) subsets such that a^^\A) is analytic in |Im -z| < 1 for 
any A^T^. Note that (jiHl) is well-defined as as-f±{A) £ J^^es C^A £ T). 

Let ^ iM^)\s=o ^('^-)) and S^iA) ^ j-/^\A)l^^ 

{A G D{6^)), where D{6u)) and D{S^) are (dense) subsets where the cor- 
responding derivatives exist. Suppose that W £ D{5S} and there exists a 
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dense subset Q± of such that 

/ dt\\[Tt{6UW)),A]\\<+^ (4-2) 

JRT 

holds for any A G Q±, where R~ (R"*") stands for the set of non-positive 
(non-negative) real numbers, then one has 



5t{A) = 6UA) + ^ dt [n (SUW)^ ,a\ , CA g Di6^)ng±) . (4-3) 



Before going to the proof of Proposition 6, we discuss its implications. Remind that 
as{A) = e^^'-=i.fl^^(^'-^'-^')Me~*^-=^.«^'-(^^~^'-^'-)' for finite-degree-of-freedom 
systems, where and {r = L, R) are, respectively, the energy and the number 
of particles in each reservoir. Hence, dui{A) = iYl,r=L nPAiHr — fir]^r),A] and 



Uw) = -^ ^[H,Hr-^^rNr] = - p,^^^^^k-J^ _ =-Js, (4-4) 

r=L,R r=L,R * ° 



where Js is the entropy production operator discussed in Sec.^ Therefore, if Propo- 
sition 6 were applicable to finite-degree-of-freedom systems, one would have 

6+{A)=i Y Pr[{Hr-f^rNr),A]+i / dt[Tt{Js),A] . 
r=L,R ■^-'^ 

On the other hand, if the state is described by a density matrix p oc e~^, it satisifes 
a KMS condition: {Aa^^{B))p = {BA)p, where a§{B) = e'^' Be'^'^' . Hence, the 
density matrix of the steady state (•••)+, if it exists, is given by a MacLennan- 
Zubarev ensemble: 



1 r ^ /-o . 

- exp - Y f3j (Hj - /u,iV,) + / dsJsis) } . (4-5) 
j=l J-oo 



As pointed out in Sec. ^ the original proposal (|4-5p by MacLennan and Zubarev 
cannot be justified. Rather, KMS states with respect to ai^^ which is generated by 
(|4-3|) should be regarded as a precise definition of the MacLennan-Zubarev ensembles. 

Proof of Proposition 6: When the dynamics is L^-asymptotic abelian and Moller 
morphisms are invertible, we have shown the same conclusion.^^^'^''-* Since the 
present conditions are different from the previous ones, we give the outline of the 
proof in case of (• • • )+. 

Remind that {Aai{B))ioc = {BA)ioc {A,B G -^"es) holds for some (dense) subset 
J^^cs -^ros- On the other hand, the map 7+ defined in Lemma 3 has the inverse 7^^ 
on J^j-cs and the steady state is given by (A)^ = {'y^{A))ioc- Hence, one has 

(^7;'h(7+(^)))>+ = (7+(^)^i(7+(5)))/oc = (7+(i?)7+(^)))zoc = {BA)+ , (4-6) 
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for any two elements of the dense set {74|7_|_(y4) G J^^cs} C J^- Namely, the steady 
state is a KMS state with respect to ai'^^ = "f^^as"f+- This proves the first half. 

Now we consider the generator. Let 7t(A) = rf^"^(rt(A)), then, for any A G JF, 
7t(A) is differentiable and ^7t(^) = ji'yt {[T-t{W),A]) , which leads to 

^7,-V,7t(^) = J [r-t {as{W) - W) ,^r'^sltiA)] 

or 



7,-V,7t(.4) = as{A) dt' [r^t, {as{W) - W) , 7,7V,7t'(^)] 



(4-7) 



where we have used 7tr_t = , t^iCTs = agT^t , 7^ ^T^t — '^-t- Since W S D{6u. 
is assumed, for any A £ D{6u)), one has 



d 



s!^HA)^-^^'a.MA) 







= 6UA) + 1 dt' We (L{W)) , a\ . (4-8) 
s=o n-J-t ^ ^ ' ^ 

It can be shown that is indeed the genrator of \^t^(yslt\s^'R.?'^^ On the other 
hand, since Lemma 3 implies 



hm ||7-V,7n(^) - cyt{A)\\ = , (4-9) 

n^oo 

for any A £ J^, one has S'^{A) = lim„^oo J^J'^ (^), if the limit exists*^. Because of 
the assumption (|4-2|) . this is indeed the case for A G D{5i^) n Q+ and we obtain the 
desired result (Q.E.D.) 

Originally, Kubo'^'^^ used the KMS condition to show the fluctuation-dissipation 
relations. Since we have a KMS characterization of NESS, it is interesting to explore 
a relation between response and fluctuation at NESS following Ref. 30). 

Suppose that a spatially uniform electric field of strength E(t) = ^e^^^ E{ijj) 
is applied to a finite domain, then, as shown e.g., by Gavish, Imry and Yurke,^^^ its 
effect is decribed by a perturbation Hamiltonian: dt'E{t')i (/ G !F: spatially 
averaged current) and we have 

/+00 J 
^ ^e-*G(a;)^M + 0(i^(t)2) , (4-10) 

where (• • • )+,_e;(j) is the state perturbed by E{t) and the frequency-dependent con- 
ductance G{ijj) is a distribution: 

1 r° 



G{iv) = ^^(^^7—5) / dte^^H[rt{I),I])+ ■ (4-11) 



*' One can show that the hmit 14-91 1 is uniform on a finite closed s-intervaL Then, the generator 
of is the graph limit of {5ir'}n>o (cf. Theorem 3.1.28 of Ref. 27)). Also (5i"' is proved to 
have the same domain; D((5i"') = D{S^). Then, if (Si"'(yl) converges, the limit is S^{A). ' 
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When E{uj) and the Fourier transform of ([r^ (/),/])+ have appropriate smoothness 
and integrabihty, the formal calculations can be justified as shown by Ruelle.^^ Note 
that the zero-frequency limit ReG(0+) of (|4-11|) corresponds to a dc-conductance, 
but is not necessarily agrees with the differential conductance as shown in 

Ref. 25) because the former comes from a local perturbation but the latter from a 
nonlocal perturbation. 

Now we introduce a gauge transformation : ^ J^, which is formally ex- 
pressed as gif{A) = e''^^r ^r+c'^ c)ip j^^-i(Y^^ Nr+c'^ c)ip ^ Clearly, (7^ is a linear map which 
preserves product and conjugation, and its action to the generators is g<p{c) = e~*'^c, 
g(fi{cLr{f)) = e~'^^ar{f) (r = L,R). Then, we observe that, if the interaction W is 
gauge-invariant: g^piW) = W, the maps Tt, ai'^^ (cf. Proposition 6) and g^p commute 
with each other. And there exists a dense subset Tg^j^-, where, for any A G J^g^ri 

Tz{A), ai'^\A) and gz{A) are analytic on the whole complex plane of z (cf. the 
arguments of Proposition 2.5.22 of Ref. 27)). Then, we have 

Corollary 7: Let (i) / £ T^^^ n DQlf^) n D(bf^) be gauge-invariant: 
g^{i) = i, (ii) W G DiSf^) n D{5P) and 

/+00 
dt\{Tt{6i)6i) + \ < +00 , (4-12) 
-00 

dt\\[Tt{6l^\w),6I]\\ < +00 iX = N,E), (4-13) 

J —00 

where 61 = I — (/)+ stands for the fluctuation of I and 6^^^ (A = N,E) 
are the infinitesimal generators of (A = N, E) formally expressed as 

(i.e., s'^^^A) = ■j^^^J^'^ {A)\s=q). Consider the Fourier transform Si{oj) of 
the symmetrized correlation function 

1 r+00 

Siiio) = - / dte''^'{{Tt{Si) 61 + 61 Tt{6i)})+ , (4-14) 
^ J -00 

then 

(e^^ - l)Si{uj) - {e^^ + l)Re(;iu;G(u;)) 



2 



00 



r dv[n{6m^^ {i6f\6i) - I ^^[rA6f\w)),6i] 



2 ^ ^ dte-y^ d.(r,(W)4+)(z4^)(5/)-/^^|[r.(5r)(H^)),5/])) 



(4-15) 



where ^s^'* is the map ^1^^ = o"i^V_^;j^5f^g, (3 = {Pl + /3ij)/2 the average 
inverse temperature, Ap = Pl — Pr the difference of inverse temperatures, 
Z\H = PlIJ-l — PnfJ'R the affinity difference and H = (PLfJ-L + PRtJ-R)/'^ 
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the average affinity. Note that the infinitesimal generator 6'^ defined by 
6+{A) = £^i'^\A)\s=o {A E D{6p) is of order of Ap and Z\H. 

From the arguments of the previous section, we have 6^^\w) = h{jR — Jl) and 

5^^'' iyV) = fi{J^ — Jl), and, thus, the terms involving t'-integrals are correlation 

functions among three current operators. On the contrary, the terms i6^^\6i) 
(A = N, E) depend on the interaction W and, in general, do not have simple physical 
meaning. For the model of an AB ring with a dot, when w = 0, (3l = (^R = and 
the left-hand side of H4-15() is absolutely integrable, the response and correlation 
functions with respect to the average current: / = — e(J/j — Jl)/2 satisfy 

|-e-^-* {(e^'^ - l)Si{uj) - (e^^ + l)Re(;;^G(cu))} 

-%r-(-(")^-'"'(l;-;l/l*'i-<")."i)). 

= ^{|r(T,(W)»')+ + | y°^<i('(T,(W)|n-(H),M])^}+0(/lM=) 

where = fJ-L — fJ-R is the chemical potential difference. Thus, the imperfection 
of the fluctuation dissipation relation is equal to a sum of two correlation functions, 
the one between the current and the interaction W and the other among three 
currents. We believe that this relation is a first step towards a quantum analog 
to the equality between the violation of fluctuation-dissipation relation and energy 
dissipation obtained for certain classical systems by Harada-Sasa^^^ and Teramoto- 
Sasa^^^ or to a nonequlibrium extension of fluctuation-dissipation relation derived 
for classical Langevin systems by Speck and Seifert.^^^ 

Proof of Corollary 7: Integrability H4-12() guarantees the existence of Si{uj) and 
hujG{uj). By integrating an analytic function e'^^^{5lTt{5I))+ on a rectangle 
[-Ti.Ta] U [T2,T2 + ih(3] U [Ts + ih(3,-Ti + ih(3] U [-Ti + ih(3,-Ti] in the com- 
plex f-plane with [zi,2;2] a segment starting from zi and terminating at Z2, and by 
taking the limit of Ti,T2 — > -|-oo, we obtain 

/+00 _ r+oo 

dte'-\5iTt{5i))+ = e-^P^ / dte'^\T_,,^0i)Tt{5i))+ . (4-16) 
-oo J — oo 

Note that the intergals on the segments \T2,T2 + ih(3] and [— Ti -|- ihp, —Ti] vanish 
for Ti,T2 +00 because of the mixing property given by Proposition 2. On the 
other hand. Proposition 6 and the gauge invariance of /: gi^{I) = I imply 

{T_,r,^{6i)n{51))+ = {T_,r,-p{g.,^{5i))Tt{5i))+ = {Tt{Si)at\T^,^f^{g,^{6im+ 

= {rt{6mt^\5i))+ = {ni5i)6i)+ + ^ C dv{ni5i)ci:\6li5i)))+ , (4-17) 

Jo 

where 61 E is shown later. It is easy to see ^i'^^ = 7^"^^i^^7+ where 

^^^=^sr%Ms,^^dd'\A)=^^^l J^^^^^^^ Thus, 
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d'^\5i) = 7+\^Tl/2&i^s/2(^+(^^))))- Since ^^^(5/) has the same structure as 
ai~^\A), one can show 61 £ D(d^) and 

= ^{^f\si) + I fjt[n{5f\w)),6i]} 

+ ^ J°^<it[r,(4^)(Ty)),5/]} , (4-18) 

from the conditions (i)-(ii) and ()4-13|) as in the proof of Proposition 6. Combining 
(ITOl) . (jTOl) . and using Si{uj) - Re{hujG{iu)) = f^dte''^\diTt{di))+ and 

Si{uj) + Re(?ttjG(a;)) = dte^'^\Tt{5I)6I)+, we obtain the desired result. (Q.E.D) 

§5. Conclusions 

As pointed out by Ruelle^^ and clearly seen from Proposition 1, nonequlibrium 
steady states investigated so far are constructed through the scattering approach. 
In this sense, the present approach can be regarded as an extension of Landauer- 
Biittiker's approach^^^ to electronic transports in mesoscopic systems. And there 
exist a class of mesoscopic systems to which the present formalism is applicable, 
such as the Ahoronov-Bohm ring with a quantum dot. Since we have a full char- 
acterization of NESS for non-interacting systems, one may develop approximations 
such as the mean-field approximation.^^) These aspects will be discussed elsewhere. 

Before closing this article, we look through a relation between the dynamical 
reversibility and irreversible evolution towards a steady state in the sense of Propo- 
sition 1.^^) We assume that the dynamics is symmetric with respect to a time 
reversal operation, namely, there exists an antilinear operation A — > l{A) such that 
l{AB) =l{A)l{B), i{aA + B) = a* l{A) + l{B) , i'^{A) = A and L{Tt{L{A))) = r_t(A) 
{A, B £ J^, a G C). The time reversal operation on a state (• • • ) is, then, defined by 
(A)^^ = {i{A'^)). For the present model, one can choose i{ar{f)) = f dkf{—k)akr, 
l(c) = c. Then, when i{W) = W (when ip = for the model of an Ahoronov-Bohm 
ring with a dot), the system has a time-reversal symmetry. Under the assumption 
that the initial state is time-reversal symmetric {A)'J-^ = {A)ioc, let us carry out the 
following thought experiment: (i) Let the system autonomously evolve up to t = io- 
(ii) The time reversal operation is performed at t = to- And (iii) let the system 
autonomously evolve once again. Just before the time reversal operation, the system 
is in the state {A)tQ- = {TtQ{A))ioc and, just after the time reversal operation, the 
state becomes 

= (A)™ = {nMA^)))ioc = {iir.toiAnioc = {T-t,{A))loc , 

which evolves further as {A)t = {Tt-to (^))/oc- Then, the system comes back at t = 2to 
to the state just before the time reversal operation, as expected. Note that, if to > 
is large enough, the state {A)tQ- just before the time reversal operation is close to 
the steady state (^)+, while the state {A)to+ just after the time reversal operation 
is close to another steady state {A)-. In other words, the time reversal operation 
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discontinuously derives the system from a state close to (•••)+ to the one close to 
(•••)_.. On the other hand, by the 'natm'al' evolution r^, the system changes towards 
the steady state {■ ■ ■)-^. Hence, the dynamical reversibility is fully consistent with 
irreversible state evolution in the sense of Proposition 1.^^^ 

It is interesting to revisit Loschmidt's criticism to the work of Boltzmann.^^) 
Although the dynamical reversibility is consistent with the irreversible state evolu- 
tion, Loschmidt's criticism can be applied to the relative entropy S{pioc\pt)- Indeed, 
in the above thought experiment, let {Js)t be the entropy production at time t, 
then, when t is slightly larger than tQ, {Js)t is close to {Js)-{< 0*^) and should it- 
self be negative. Thus, because of ()3-22|) . there is a period when the relative entropy 
S{pioc\pt) decreases. Contrary to Boltzmann's reply, '^^^ these states are typical in the 
sense that they evolve towards the steady state (•••)+ for t oo. In other words, 
Loschmidt's criticism does not deny the consistency of irreversible phenomena with 
dynamical reversibility and it just shows that the relative entropy is not an appro- 
priate thermodynamic entropy for general cases. Another criticism to Boltzmann's 
work by Zermelo^^^ is not applicable to the present case since the recurrence time is 
infinitely long as a result of the infinite extension of the system. 

From the point of view of the second law of thermodynamics, one may be sat- 
isfied with all these features, particularly the properties of the entropy production 
discussed in Sec. 3. However, one should remind that the 'correct' form of the en- 
tropy production is obtained because we start from a local equilibrium state where 
each subsystem is in a canonical state. As the canonical states are very outcome of 
the second law, the present results do not give a dynamical proof of the second law, 
but show that, once one starts from canonical ensembles or their combinations, the 
dynamics derives the system consistently with the second law of thermodynamics. 
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